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Abstract 

Dijet cross sections as functions of several jet observables are measured in photoproduction 
using the HI detector at HERA. The data sample comprises e'^p data with an integrated 
luminosity of 34.9 pb^^. Jets are selected using the inclusive k± algorithm with a minimum 
transverse energy of 25 GeV for the leading jet. The phase space covers longitudinal 
proton momentum fraction Xp and photon longitudinal momentum fraction in the ranges 
0.05 < Xp < 0.6 and 0.1 < x-^ < 1. The predictions of next-to-leading order perturbative 
QCD, including recent photon and proton parton densities, are found to be compatible with 
the data in a wide kinematical range. 
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1 Introduction 



In QCD (Quantum Chromo Dynamics) the photoproduction of jets with high transverse energy 
is described by the hard interaction of real photons with quarks and gluons inside the proton. In- 
teractions with two outgoing partons of large transverse momentum are due to direct processes, 
such as 7g — > gq (QCD-Compton effect) and 7(7 — > qq (photon-gluon fusion) and resolved 
processes where the photon first splits into a quark pair (or higher multiplicity fluctuation) and 
one of the resulting partons subsequently scatters off a parton in the proton. The calculation of 
the latter processes can be approximated by ascribing parton densities to the photon, which also 
include the inherently non-perturbative aspects of the photon structure. 

In analogy to the proton case parton densities of the photon depend on a factorization scale 
fij and on x^, the longitudinal momentum fraction of the photon taken by the interacting par- 
ton. The limiting case of direct interactions is given by = 1. At HERA these photopro- 
duction reactions can be investigated in inelastic electron (positron) proton reactions at very 
small squared four-momentum transfers Q^. Starting from the first investigation of this kind at 
HERA [jT]] the comparison of the predictions of QCD with the results has been a central topic of 
interest [0-^. These investigations are particularly interesting, because previous measurements 
of high transverse energy jet production in ep and pp scattering were not fully described by 
QCD calculations [g-^]. 

High transverse energy jets provide a natural hard scale for perturbative QCD calculations. 
Such calculations have been performed for direct and resolved processes in leading (LO) and 
next-to-leading (NLO) order. The measurement of jet cross sections at high transverse energy 
presented in this paper can therefore be used to test the current predictions of NLO perturbative 
QCD and the parameterizations of photon and proton parton densities at large scales with a 
precision of typically 10%. Photon quark densities have been determined in experiments at 
e"'"e~ -colliders which investigate the photon structure function F2, where x-y values up to 
0.8 and scales up to 500 GeV^ have been reached. In comparison the analysis presented here 
extends the range up to 1 at scales between 600 and 6000 GeV^, where the quark density 
parameterizations of the photon are presently not well constrained by measurements. In contrast 
to the F2 measurements, the photoproduction of jets is directly sensitive to the gluon density 
of the photon, which is poorly known to date. Furthermore our data are sensitive to the parton 
densities of the proton at fractional momentum values Xp up to 0.6. In this kinematical regime, 
the quark densities are well known from deeply inelastic scattering data, while the gluon density 
has uncertainties of the order 10 to 50% [|^]. Photoproduction data can thus be used to constrain 
the parton density functions in regions where only few measurements are presently available. 
However, detailed parton densities can not be extracted from these data alone. 

This paper is based on an e^p data sample collected with the HI detector in the years 1995- 
1997 and corresponds to an integrated luminosity of 34.9 pb^^. It presents dijet cross sections 
as a function of jet observables with mean jet transverse energies 20 < Ex < 80 GeV, observed 
x^ values 0.1 < x^ < 1 and values of Xp ranging from 0.05 to 0.6. The two jets considered in 
the investigated process ep e jet jet X are defined as the two jets with the highest transverse 
energyQ 

'a similar analysis has recently been made available |jl0|]. 
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2 Jets in Photoproduction 



2.1 Cross sections and observables 



The cross section for the photoproduction of hard jets in electron-proton collisions, a^p, can be 
calculated from the photon-proton scattering result, a^p, using the factorization ansatz 



Here the usual variable y of deeply inelastic scattering is interpreted as the longitudinal momen- 
tum fraction of the incoming electron taken by the photon and g is the photon flux calculated 
in the Weizsacker- Williams approximation lO-O]. The hadronic photon-proton jet cross sec- 



tion is obtained as the convolution of the partonic cross sections with the parton momentum 
distributions of the proton /j/p and the photon fjj^ . As outlined in the introduction it is divided 
into a sum of two components, the direct part a'^'p', where the photon directly interacts with a 
parton of the proton and the resolved part c^^p°''^^ where one of the partons inside the photon 
interact with a parton of the proton. This distinction is unambiguously defined in leading order 
only and depends on the photon factorization scale /x^. The two components can be expressed 
as: 



^direct 



a. 



resolved 
IP 



The squared centre-of-mass energy of the hard subprocess is s = XpX-^ys, where ^/s is the total 
centre-of-mass energy in the ep-system, i.e. 300 GeV for this analysis. The proton factorization 
scale is jjip and the renormalization scale is /i^. The partonic cross sections a can be expanded 
as a perturbative series in powers of as and have been calculated up to the next-to-leading order 
in QCD[lTf^]. 

The total cross sections on the left hand side of equation [l] are obtained by integrating over y, 
Xp and x^. The partonic cross sections and (Jij contain a further integration over an internal 
degree of freedom, e.g. cos 6* , the scattering angle in the centre-of-mass system of the partonic 
two body reaction, or the transverse energy. More detailed information on the reaction dynamics 
is obtained by measuring differential cross sections in these kinematical observables. In order 
to avoid singularities in the partonic cross sections a minimum cut in 6* or in the transverse 
energy of the outgoing partons has to be applied. 

The two scaled longitudinal parton momenta x^ and Xp are calculated from the jets produced 
in the hard subprocess, using the definition 

X, = -^{ET,ie-^^+ET.2e"^') (4) 

Xp = -^{ET,ie^^ + ET,2e'i') (5) 



Here Ex^i and Et.2 are the transverse energies of the two jets of the hard subprocess, rji and 772 
are their pseudorapidities in the laboratory frame (r] = — ln(tan6'/2)) and E^. and Ep are the 
energies of the electron and proton beams[[ These relations are used as definitions of observ- 
ables in all orders and are easily derived for 2 ^ 2 processes, where the transverse energies of 
the jets are equal. The pseudorapidities of the jets are related to 9* via 

cos^* = I tanh((r/i - V2)/2)\ . (6) 

In principle one could measure the dependence of the fourfold differential cross section 
cos9* on all four variables. This, however, would require a much larger data 
set than presently available. Therefore in this paper more inclusive quantities are presented. The 
distribution of the invariant mass of the two jets with the highest transverse energies, Mjj, the 
mean transverse energy of the two leading jets ET,mean and the transverse energy distribution 
of the highest transverse energy jet, Ex^max, are studied. The cross section differential in the 
average value of the pseudorapidities = (r/i + ?72)/2 is particularly sensitive to parton density 
functions. It is thus presented for different photon-proton centre-of-mass energies (y regions) 
and different scales (Ex^max regions), cf. equations ^ and ^. 

Differential cross sections in and Xp are measured in different scale regions {E^^max re- 
gions) and for different x^ or Xp cut-off values. The angle 6* is sensitive to the dynamics of jet 
production and the corresponding differential cross section is therefore evaluated for different 

regions for all Mjj and in addition with a cut in Mjj. The 0086** distribution could be in- 
fluenced by the production of W or bosons, whose hadronic decays have a different angular 
distribution from that expected for QCD dijet production. Using the EPVEC Monte Carlo gen- 
erator [[TSl], the contribution of W bosons is estimated to be 5-6 events. The contribution is 
expected to be negligible. The background from these processes is therefore not considered in 
the following. 

In the present analysis jets are defined using the inclusive algorithm as proposed in 
[]r9|,pD|]. The application of this algorithm has become standard in jet analyses at HERA pT|]. 
It utilizes a definition of jets in which not all particles are assigned to hard jets. Here it is 
applied in the laboratory frame with the separation parameter set to 1 and using an E^ weighted 
recombination scheme in which the jets are considered massless. 



2.2 QCD Predictions and Models 

To simulate the direct and resolved photoproduction of jets, the PYTHIA 5.7 p2|] and HERWIG 



5.9 [ P3| ] event generators were used followed by a full detector simulation [ P4{ | of all Monte 
Carlo events. Both programs contain the Born level QCD hard scattering matrix elements, 
regulated by a minimum cut-off in transverse momentum. Leading logarithmic parton showers 
are used to represent higher order QCD radiation. GRV-LO parton density functions 

(pdfs) for the proton and photon were chosen. The Lund String model is applied in PYTHIA to 
hadronize the outgoing partons, while in HERWIG the cluster hadronization approach is used. 

^The coordinate system is centered at the nominal interaction point with the positive z direction along the 
incident proton beam. The polar angle is defined with respect to the positive z axis. 
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Multiple interactions between the proton and the resolved photon are dealt with in PYTHIA by 
adding additional interactions between spectator partons within the same event. These processes 
are calculated by extending the perturbative parton-parton scattering to a low Et cut-off. 

In HERWIG multiple interactions are included by producing in a fraction P' of the resolved 
events so called soft underlying events. These interactions are parameterized using experimental 
results of soft hadron-hadron scattering. The effect of multiple interactions is tested by com- 
paring, in the data and in the HERWIG calculations, the energy flow distributions around the 
jet axis with and without a fraction P' of events containing the soft underlying event. For P' 
~30-35% these distributions are found to be well described for all regions of Xj. The differ- 
ence of the calculated HERWIG cross sections with and without 35% of soft underlying events 
is below 10% for between 0.3-0.8 and 10-20% for x^ < 0.3. For Xj > 0.8 the difference is 
negligible. PYTHIA is also able to describe these distributions. 

The goal of this analysis is the comparison of the measured cross sections to perturbative 
QCD calculations at the parton level. The LO and NLO dijet cross sections were computed 



using a program based on the subtraction method [ p^pTQ for the analytic cancellation of infrared 
singularities. In calculating LO and NLO cross sections a 2-loop as was taken with 5 active 
flavours. Aqcd was set to 0.226 GeV (as{Mz) = 0.118), which is the value used in the proton 
parton density functions. CTEQ5M parton density functions were chosen for the proton 
whereas MRST99 p9|] parton density functions were selected to test the dependence of the 



NLO cross sections on the proton pdfs. For the photon we choose GRV-HO [ pOQ as a main 
setting and the parameterization of AFG-HO [ PT| ] to study the dependence of the results on the 
choice of the photon pdfs. The renormalization scale fir and the factorization scales fip and ji^ 
were, event by event, set to the sum of the transverse energies of the outgoing partons divided 
by two. The QCD program allows the variation of this common scale. It was varied from 
0.5 to 2 times the default scale to estimate the scale uncertainty in the NLO calculation. This 
uncertainty turned out to vary between ±10 and ±20% in the measured kinematic range. 

In addition the data are compared to the predictions of NLO QCD corrected for hadroniza- 
tion effects, which are defined as the ratio of the cross sections with jets reconstructed from 
hadrons and from partons before hadronization. The hadronization effects are calculated with 
PYTHIA and HERWIG and the mean value of the two predictions is used for corrections. Here 
the difference between the two Monte Carlo models is in general very small and at maximum 
10%. The jets built out of partons are found to be very well correlated with the jets built out of 
hadrons. 



3 Experimental Technique 
3.1 HI Detector 

The HI detector is described in detail in Only those components relevant to the present 



analysis are briefly described here. The Liquid Argon (LAr ) [ |34| ] and SpaCal [ p5| | calorimeters 
were used to trigger events, to reconstruct the hadronic energy of the final state and to select 
photoproduction events by eliminating events with an identified scattered positron. The LAr 



calorimeter covers the polar angle range A° < 6 < 154° with full azimuthal acceptance. The 
jet energy calibration agrees at the 2% level with the Monte Carlo simulation as determined 
by the transverse energy balance between jet and electron for deeply inelastic scattering events 
and by the transverse energy balance between the two jets for the photoproduction sample in 
different kinematic regions. The angular region 153° < 9 < 177.8° is covered by the SpaCal, 
a lead/scintillating-fibre calorimeter. It has a hadronic energy scale uncertainty of 8%. The 
central tracking detector (CJC) was used to reconstruct the interaction vertex and to supplement 
the measurement of hadronic energy flow. The CJC consists of two concentric cylindrical drift 
chambers, coaxial with the beam-line, with a polar angle coverage of 15° < 9 < 165°. The 
entire CJC is immersed in a 1.15 T magnetic field. The luminosity determination is based on 
the measurement of the ep ep'j Bethe-Heitler process, where the positron and photon are 
detected in calorimeters located downstream of the interaction point in e-beam direction. 



3.2 Event Selection 

The data sample was collected at HERA with the HI detector in the years 1995-97, when 
protons of 820 GeV energy collided with positrons of 27.6 GeV energy resulting in a centre- 
of-mass energy of 300 GeV. The events were triggered on the basis of high transverse energy 
deposits in the LAr calorimeter. The trigger efficiencies were above 94% for the event sample 
described in this analysis. Energy deposits in the calorimeters and tracks in the CJC were com- 
bined in a manner that avoids double counting to reconstruct the hadronic energy of events p^ ] . 

It was required that an event vertex was reconstructed within 35 cm of the nominal z po- 
sition of the vertex. The most significant background in the data sample arises from neutral 
current deeply inelastic scattering events, and was suppressed by removing events with an elec- 
tron identified in the LAr calorimeter or SpaCal and by requiring y < 0.9, with y reconstructed 
using hadronic variables [B7|]. This reduces the background to less than 1% for the total sam- 



ple. In the region with the highest y at low f] the remaining background was calculated to be 
about 5% based on deeply inelastic scattering dijet data and the ARIADNE [ pS] ] Monte Carlo 
interfaced with DJANGO [|39[]. It was subtracted statistically. After applying a cut on the 
missing transverse energy ET^miss < 20 GeV the remaining charged current (ep uX) and 
non-ep scattering background was found to be negligible. Events induced by cosmic rays were 
removed. 

Asymmetric cuts on the Et of the two jets with the highest transverse energies are applied to 
avoid regions of phase space affected by uncertainties in the NLO calculation [ [l4| ] . On the other 
hand a highly asymmetric cut causes large NLO corrections and a pronounced dependence on 
the choice of scale. The jet selection criteria therefore required an Et of the highest transverse 
energy jet Ex^max > 25 GeV, and the transverse energy of the second highest transverse energy 
jet Et ^second > 15 GeV. When the cut on E^. second is varied between ±5 GeV, the ratio of the 
measured cross sections to the theoretical prediction varies by up to 10% for < 0.8 and by 
up to 3% for > 0.8. 

The pseudorapidity of each jet r/j was restricted to —0.5 < rji < 2.5. All jets are thus well 
contained in the LAr calorimeter. The measured kinematic region was restricted to 0.1 < y < 
0.9 and Q"^ < \ GeV^, as given by the acceptance for electrons in the LAr and Spacal. The 
kinematic range of the measured dijet cross sections is summarized in Table [1]. Applying these 
cuts the total number of events measured was 5265. 
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< 1 GeV" 
0.1 < y < 0.9 
> 25 GeV 

Et, second > 15 GeV 

-0.5 <rji< 2.5 

Table 1: The definition of the phase space of the measured dijet cross sections. 



3.3 Correction of the Data for Detector Effects 

The data were corrected for detector effects such as limited resolution and inefficiencies. To 
determine these effects the HERWIG and PYTHIA Monte Carlo samples were used. Both 
programs do not describe the absolute normalization of the dijet cross sections. After scaling 
the HERWIG cross sections by 2 and the PYTHIA cross sections by 1.2 the two programs gave 
a good description of the measured jet distributions in shape and normalization. 

The bin sizes of all distributions are matched to the resolution and result in good bin ef- 
ficiency and purity. The correction was done by using the so called bin-to-bin method. The 
correction functions were calculated from the ratio of the cross sections with jets reconstructed 
from hadrons (hadron level) and from detector objects (detector level) in each bin, where each 
sample was subject to the selection criteria defined above. The correction functions of the two 
models are in good agreement and differ on average by 5% and at most by 20%. The mean 
values of the two Monte Carlo generators were thus taken for the correction. The resulting 
correction factors typically have values between 0.8 and 1.2. 

3.4 Systematic Uncertainties 

For the jet cross sections the following sources of systematic error were considered: 

• A 2% uncertainty in the LAr energy scale results in an uncertainty of typically 10%. 

• An 8% uncertainty in the hadronic Spacal energy scale results in an uncertainty of 1%. 

• In addition to the variations of the calorimeter energy scales a shift of 1% on y is consid- 
ered. This variation results in an uncertainty of 3%. 

• Half of the difference between the correction factors calculated with HERWIG and with 
PYTHIA is taken as the uncertainty in the detector correction. The resulting uncertainty 
is less than 10%. 

• The uncertainty in the trigger efficiency results in an error of ~ 3%. 

• The uncertainty in the background subtraction results in an error of ~ 2%. 

• The uncertainty in the integrated luminosity results in an overall normalization error of 
1.5%. 
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The statistical and all systematic errors are added in quadrature. The resulting total uncer- 
tainty ranges from 10 to 30%, where the systematic contribution is dominated by uncertainties 
in the calorimeter energy scales and in the correction to the hadron level. 



4 Results 

The measured cross sections for inclusive dijet production in the reaction ep — > e jet jet X are 
given as single differential cross sections in all cases. The data are corrected for detector effects 
and are presented at the level of stable hadrons for the phase space region defined in Table [l]. 
The inner error bars of the data points in the figures denote the statistical, the outer error bars 
the total uncertainty. The data are also presented in Tables All results are compared to 
next-to-leading order (NLO) QCD predictions obtained with the standard setting described in 
section ^ if not otherwise quoted. The predictions of NLO QCD corrected for hadronization 
effects NL0(1 + 5hadr) are also shown. 

In Figure |l] the dijet cross section is shown as a function of the invariant mass Mjj of the 
dijet system. The data are presented for Mjj values between 45 and 180 GeV. The measured 
cross section falls by about 3 orders of magnitude over this range. NLO QCD describes the 
measured cross sections for the whole mass range. Hadronization corrections are less than 5% 
for all bins. The calculation using LO matrix elements fails to describe the low Mjj region. This 
is partly due to the fact that the low Mjj region is populated by events which are influenced 
by the asymmetric cuts on the jet transverse energies. Events in which the second jet has a 
transverse energy below 25 GeV contribute mainly in this region. In dijet calculations they 
only appear beyond leading order. The scale uncertainties in the QCD predictions are largest at 
low Mjj values. 

A similar statement on the large scale uncertainties and the difference between data and 
the LO calculation holds for small transverse momenta. In Figure ^) the dijet cross section 
da / dET,mean is shown. Here the scale uncertainties decrease from ±20% for the first bins to 
less than ±5% for ET,mean > 30 GeV. The data are well described by the NLO calculation. 
The dijet cross section as a function of the transverse energy of the highest transverse energy jet 
ET,max is shown in Figure]^). The distribution again demonstrates that the data are described 
by NLO QCD up to the highest Er^max values within errors. The NLO scale uncertainty is 
not reduced significantly with increasing Ex^max- The cross sections differential in transverse 
energy are hardly altered by hadronization corrections which are around 5% for all bins. The 
NLO QCD calculation with hadronization corrections predicts the measured cross sections up 
to the highest masses and transverse energies, although the photon and proton pdfs have been 
extracted from quite different processes and mostly at lower scales. 

To further explore the photon and proton structure the differential cross section da/drj is 
displayed in Figure ^for two ranges of Ex,max subdivided into two y regions. While the former 
implies a variation of the scale the latter corresponds to different center of mass energies in the 
photon-proton-system. Again, good agreement between data and NLO QCD is observed taking 
into account the uncertainties in the calculations and in the data points. The predictions tend 
to lie above the data at low f], where direct interactions dominate and hadronization corrections 
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are largest. At high rj, where in contrast resolved interactions dominate and hadronization 
corrections are small, the NLO QCD predictions agree well with the measured data. 

Figures ^a) and b) show the dijet cross section da/dx^ as a function of for two different 
Xp regions. The calculations exceed the data, while remaining within the given uncertainties, 
only for > 0.85, where the largest hadronization correction occur. Using the MRST99 1-3 
proton pdfs (with a large variation of the high Xp gluon density) instead of CTEQ5M results 
in differences of less than 5% for the predicted cross section for Xp < 0.1 and up to 15% for 
Xp > 0.1. This is smaller than the scale uncertainties for Xp < 0.1 and of the same order for 
Xp > 0.1. These findings are corroborated in Figures ^) and d) where the cross section da/dxp 
is shown as a function of Xp for two different x^ regions. Even at the highest Xp the measured 
cross sections are seen to agree well with the QCD predictions, which in this part of the phase 
space attribute about 40% of the cross section to processes induced by gluons in the proton. 
The constraints on the pdfs used in the QCD calculations here come dominantly from deeply 
inelastic scattering at lower scales where the gluon fraction is smaller. The concept of universal 
pdfs in hard processes in QCD is thus observed to describe measurements with rather different 
experimental conditions. 

Figure |] displays the dijet cross sections da/dx-y as a function of x^ for two regions of 
ET,max, representing different factorization scales for the photon and proton pdfs. The data are 
compared to NLO calculations corrected for hadronization effects with two different parame- 
terizations of the photon structure. The predictions describe the data well and vary only slightly 
with the photon pdfs used. In contrast the NLO scale uncertainties produce a significant effect 
as can be inferred from Figure H, which repeats the data of Figure ^ with a comparison of the 
GRV-HO pdfs of the photon. For high values of x^ the hadronization corrections are sizeable 
and improve the agreement with the data. A more detailed comparison between data and theory 
is obtained by plotting their relative difference as shown in Figure 0. NLO predictions including 
hadronization corrections are shown for both sets of photon pdfs. At variance to the previous 
plots the error bars of the data contain only the uncorrelated systematic errors, while the corre- 
lated errors due to the uncertainty in the calorimeter energy scales are shown as a hatched band. 
Figure ^ shows that the assumed NLO scale uncertainties are the dominant source of uncertain- 
ties in the comparison of data and theory. The expectation exceeds the data only for the high x^ 
and high Ex^max regions. Within these uncertainties the picture of an universal photon structure 
is thus corroborated. 

Finally, the dijet cross section da /d cos 6* is plotted in Figures ^) and b) for x^ < 0.8 and 
x^ > 0.8 respectively. The cross section decreases with increasing cos6'* mainly because of 
the cuts in Ex- Again, the data are well described by NLO QCD for low x^, whereas at higher 
Xj the predictions overshoot the data for small values of cos^^*. These cross sections are also 
shown with a cut on the invariant mass Mjj of the dijet system in Figures ^) and d), essentially 
excluding the first bin of Figure [1|. The cut reduces the restriction of the phase space due to 
the correlation with the Et requirements and changes the shape of the distribution towards that 
expected from the QCD matrix elements. The QCD calculations reproduce this transition nicely 
in both x^ regions where resolved and direct photon induced processes contribute with different 
weights. 
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5 Conclusions 



New measurements of dijet cross sections in photoproduction at high transverse energies are 
presented for various jet kinematic observables. The measurements cover invariant dijet masses 
up to 180 GeV and transverse energies up to 80 GeV, reaching Xp and values where the 
experimental information was previously limited. In this kinematic domain non-perturbative 
effects like multiple interactions and hadronization are found to be small, which allows a direct 
comparison of NLO QCD calculations with the data to be made. The results demonstrate the 
power of perturbative QCD in predicting the measured cross sections in a wide kinematical 
range. Even though the photon pdfs have been obtained from measurements at lower scales, 
their QCD evolution correctly reproduces the data at high scales. The data do not require sig- 
nificant changes in the parameterizations of the pdfs but are certainly useful to further constrain 
the existing ones. Likewise our understanding of the proton structure in the high Xp, high scale 
region can be improved with the help of these data. A future stronger constraint requires a 
reduction of both the theoretical scale uncertainties and the systematic uncertainties in the data. 
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33 


32 



Table 2: Differential ep cross sections for dijet production as a function of the invariant dijet 
mass Mjj (upper table), as a function of Ex^mean (middle table) and as a function of E-r^max 
(lower table) with statistical and total upper and lower uncertainties (cf. Figures [T] and|^). 
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Table 3: Differential ep cross sections for dijet production as a function of rj with statistical and 
total upper and lower uncertainties (cf. Figure |]). 
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Table 4: Differential ep cross sections for dijet production as a function of x^ with statistical 
and total upper and lower uncertainties (cf. Figure 
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Table 5: Differential ep cross sections for dijet production as a function of Xp with statistical 
and total upper and lower uncertainties (cf. Figure 
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Table 6: Differential ep cross sections for dijet production as a function of x^ with statistical 
and total upper and lower uncertainties (cf. Figures |], || and 0). 
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Table 7: Differential ep cross sections for dijet production as a function of cos 6* with statistical 
and total upper and lower uncertainties (cf. Figure |8|). 
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Figure 1: Differential ep cross sections for dijet production (Q^ < IGeV^) as a function of the 
invariant dijet mass Mjj of the two highest Et jets. Here, as well as in the following figures 
unless explicitly stated otherwise, the inner error bars denote the statistical error, the outer er- 
ror bars the total uncertainties of the data. The LO predictions using CTEQ5M pdfs for the 
proton and GRV-HO pdfs for the photon are shown as a dotted line. NLO predictions with 
the same pdfs are shown as a dashed line. The full line shows the NLO predictions, includ- 
ing hadronization corrections and the grey band indicates the renormalization and factorization 
scale uncertainties of the NLO prediction. 
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Figure 2: Differential ep cross sections for dijet production (Q^ < IGeV^) as a function of 
a) ET,mean, the mean and b) ET,max, the maximum Et of the two highest Et jets. The LO 
predictions using CTEQ5M pdfs for the proton and GRV-HO pdfs for the photon are shown as 
a dotted line. NLO predictions with the same pdfs are shown as a dashed line. The full line 
shows the NLO predictions, including hadronization corrections and the grey band indicates the 
renormalization and factorization scale uncertainties of the NLO prediction. 
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Figure 3: Differential ep cross sections for dijet production (Q^ < IGeV^) as a function of the 
f] of the two highest Et jets. The regions of low y, a) and b) and high y, c) and d) are shown for 
two ranges of Ex^max- NLO predictions using CTEQ5M pdfs for the proton and GRV-HO pdfs 
for the photon are shown as a dashed line. The full line shows the NLO predictions, includ- 
ing hadronization corrections and the grey band indicates the renormalization and factorization 
scale uncertainties of the NLO prediction. 
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Figure 4: Differential ep cross sections for dijet production (Q^ < IGeV^) as a function of a) 
and b) and Xp c) and d). Figures a) and b) distinguish regions of small and large Xp and figures 
c) and d) corresponding regions in x-y. NLO predictions using CTEQ5M pdfs for the proton 
and GRV-HO pdfs for the photon are shown as a dashed line. The full line shows the NLO 
predictions, including hadronization corrections and the grey band indicates the renormalization 
and factorization scale uncertainties of the NLO prediction. 
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Figure 5: Differential ep cross sections for dijet production (Q^ < IGeV^) as a function of 
Xj for a) low ET^max and b) high Ex^max- The NLO predictions using CTEQ5M pdfs for the 
proton and GRV-HO pdfs for the photon and including hadronization corrections are shown as 
a full line. NLO predictions using AFG-HO parametrizations of the photon pdfs and including 
hadronization corrections are shown as the dashed line. 
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Figure 6: Differential ep cross sections for dijet production (Q^ < IGeV^) as a function of 
for a) low Ex^max and b) high Ex^max- NLO predictions using CTEQ5M pdfs for the proton 
and GRV-HO pdfs for the photon are shown as a dashed line. The full line shows the NLO 
predictions, including hadronization corrections and the grey band indicates the renormalization 
and factorization scale uncertainties of the NLO prediction. 
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Figure 7: The dependence of the relative difference of the measured dijet cross sec- 
tions (Q^ < IGeV^) from the NLO prediction, with hadronization corrections applied using 
CTEQ5M pdfs for the proton and GRV-HO pdfs for the photon (here axheory)- The symbol 
a stands for da/dx^. Shown is the relative difference of the data (points) and the NLO pre- 
dictions using the AFG-HO pdf (dashed line) with hadronization corrections applied. Figures 
a) and b) show the relative difference for the lower Ex^max and higher Ex^max regions respec- 
tively. The inner error bars denote the statistical error, the outer error bars denote all statistical 
and uncorrected systematic errors of the data added in quadrature. The correlated systematic 
errors are shown in the middle plots as a shaded band. The grey band (lower plots) shows the 
renormalization and factorization scale uncertainties of this NLO prediction. 
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Figure 8: Differential ep cross sections for dijet production (Q^ < IGeV^) as a function of 
cos6'* distinguished for small Xj a) and c) and large Xj b) and d). Figures c) and d) show the 
cross sections for large invariant masses of the dijet system. NLO predictions using CTEQ5M 
pdfs for the proton and GRV-HO pdfs for the photon are shown as a dashed line. The full line 
shows the NLO predictions, including hadronization corrections and the grey band indicates the 
renormalization and factorization scale uncertainties of the NLO prediction. 
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